Search for: All records

We report the discovery of two new Galactic accreting compact objects consistent with the respective positions of the unassociated Fermi-LATgamma-ray sources 4FGL J0639.1-8009 and 4FGL J1824.2+1231. A combination of new and archival X-ray data from Chandra, XMM-Newton, Swift/XRT, and eROSITA reveals a variable X-ray source in each gamma-ray error ellipse. Both candidate counterparts show power-law spectra with photon indices Gamma ∼ 1.7 to 1.9. Optical follow-up photometry and spectroscopy show rapid high-amplitude variability unrelated to orbital motion and persistent accretion disk spectra for both objects. We demonstrate that the properties of these X-ray/optical sources are at odds with the known phenomenology of accreting white dwarfs, but are consistent with the observed properties of the subluminous disk state of transitional millisecond pulsars. This brings the census of confirmed or candidate transitional millisecond pulsars in the Galactic field to nine. We show this potentially represents <~ 10% of the total population of transitional millisecond pulsars within 8 kpc. 

ABSTRACT While nova eruptions produce some of the most common and dramatic dust formation episodes among astrophysical transients, the demographics of dust-forming novae remain poorly understood. Here, we present a statistical study of dust formation in 40 novae with high-quality optical/IR light curves, quantitatively distinguishing dust-forming from non-dust-forming novae while exploring the properties of the dust events. We find that 50–70 per cent of novae produce dust, significantly higher than previous estimates. Dust-forming novae can be separated from those that do not show dust formation by using the largest redward ($V-K$) colour change from peak visible brightness; ($V-J$) or ($V-H$) offer useful but less sensitive constraints. This makes optical+IR photometry a powerful tool to quantify dust formation in novae. We find that novae detected in GeV $$\gamma$$-rays by Fermi-LAT appear to form dust more often than novae not detected by Fermi, implying a possible connection between $$\gamma$$-ray-producing shocks and dust production. We also find that novae that evolve very quickly ($$t_2 < 10$$ d) are much less likely to form dust, in agreement with previous findings. We confirm a correlation between $$t_2$$ and the time of the onset of dust formation (which occurs $$\sim$$1 week–3 months after maximum light), but conclude that it is primarily an observational artefact driven by dust formation determining when a nova drops 2 mag below peak. The significant fraction of novae that form dust make them ideal laboratories in our Galactic backyard to tackle the puzzle of dust formation around explosive transients. 

Abstract We present a systematic study of the BVRI colours of novae over the course of their eruptions. Where possible, interstellar reddening was measured using the equivalent widths of Diffuse Interstellar Bands (DIBs). Some novae lack spectra with sufficient resolution and signal-to-noise ratios; therefore, we supplement as necessary with 3D and 2D dust maps. Utilising only novae with DIB- or 3D-map-based E(B − V), we find an average intrinsic (B − V)0 colour of novae at V-band light curve peak of 0.20 with a standard deviation of 0.31, based on 25 novae. When the light curve has declined by 2 magnitudes (t2), we find an average (B − V)0 = −0.03 with a standard deviation of 0.19. These average colours are consistent with previous findings, although the spreads are larger than previously found due to more accurate reddening estimates. We also examined the intrinsic (R − I)0 and (V − R)0 colours across our sample. These colours behave similarly to (B − V)0, except that the (V − R)0 colour gets redder after peak, likely due to the contributions of emission line flux. We searched for correlations between nova colours and t2, peak V-band absolute magnitude, and GeV γ-ray luminosity, but find no statistically significant correlations. Nova colours can therefore be used as standard “crayons” to estimate interstellar reddening from photometry alone, with 0.2–0.3 mag uncertainty. We present a novel Bayesian strategy for estimating distances to Galactic novae based on these E(B − V) measurements, independent of assumptions about luminosity, built using 3D dust maps and a stellar mass model of the Milky Way. 

ABSTRACT V745 Sco is a Galactic symbiotic recurrent nova with nova eruptions in 1937, 1989, and 2014. We study the behaviour of V745 Sco at radio wavelengths (0.6–37 GHz), covering both its 1989 and 2014 eruptions and informed by optical, X-ray, and $$\gamma$$-ray data. The radio light curves are synchrotron-dominated. Surprisingly, compared to expectations for synchrotron emission from explosive transients such as radio supernovae, the light curves spanning 0.6–37 GHz all peak around the same time ($$\sim$$18–26 d after eruption) and with similar flux densities (5–9 mJy). We model the synchrotron light curves as interaction of the nova ejecta with the red giant wind, but find that simple spherically symmetric models with wind-like circumstellar material (CSM) cannot explain the radio light curve. Instead, we conclude that the shock suddenly breaks out of a dense CSM absorbing screen around 20 d after eruption, and then expands into a relatively low-density wind ($$\dot{M}_{out} \approx 10^{-9}\!-\!10^{-8}$$ M$$_{\odot }$$ yr$$^{-1}$$ for $$v_w = 10$$ km s$$^{-1}$$) out to $$\sim$$1 yr post-eruption. The dense, close-in CSM may be an equatorial density enhancement or a more spherical red giant wind with $$\dot{M}_{in} \approx [5\!-\!10] \times 10^{-7}$$ M$$_{\odot }$$ yr$$^{-1}$$, truncated beyond several $$\times 10^{14}$$ cm. The outer lower-density CSM would not be visible in typical radio observations of Type Ia supernovae: V745 Sco cannot be ruled out as a Type Ia progenitor based on CSM constraints alone. Complementary constraints from the free–free radio optical depth and the synchrotron luminosity imply the shock is efficient at accelerating relativistic electrons and amplifying magnetic fields, with $$\epsilon _e$$ and $$\epsilon _B \approx 0.01\!-\!0.1$$. 

Abstract We analyse stellar streams in action-angle coordinates combined with recent local direct acceleration measurements to provide joint constraints on the potential of our galaxy. Our stream analysis uses the Kullback–Leibler divergence with a likelihood analysis based on the two-point correlation function. We provide joint constraints from pulsar accelerations and stellar streams for local and global parameters that describe the potential of the Milky Way (MW). Our goal is to build an “acceleration ladder,” where direct acceleration measurements that are currently limited in dynamic range are combined with indirect techniques that can access a much larger volume of the MW. To constrain the MW potential with stellar streams, we consider the Palomar 5, Orphan, Nyx, Helmi, and GD1 streams. Of the potential models that we have considered here, the preferred potential for the streams is a two-component Staeckel potential. We also compare the vertical accelerations from stellar streams and pulsar timing, defining the function $f\left(z\right)={\alpha }_{1\mathrm{pulsar}}z-\frac{\partial \mathrm{\Phi }}{\partial z}$, where Φ is the MW potential determined from stellar streams andα_{1 pulsar}zis the vertical acceleration determined from pulsar timing observations. Our analysis indicates that the Oort limit determined from streams is consistently (regardless of the choice of potential) lower than that determined from pulsar timing observations. The calibration we have derived here may be used to correct the estimate of the acceleration from stellar streams. 

ABSTRACT We present a model for the formation of the Magellanic Stream (MS) due to ram pressure stripping. We model the history of the Small and Large Magellanic Clouds in the recent cosmological past in a static Milky Way (MW) potential with diffuse halo gas, using observationally motivated orbits for the Magellanic Clouds derived from HST proper motions within the potential of the MW. This model is able to reproduce the trailing arm but does not reproduce the leading arm feature, which is common for models of the stream formation that include ram pressure stripping effects. While our model does not outperform other models in terms of matching the observable quantities in the MS, it is close enough for our ultimate goal – using the MS to estimate the MW mass. By analysing our grid of models, we find that there is a direct correlation between the observed stream length in our simulations and the mass of the MW. For the observed MS length, the inferred MW mass is 1.5 ± 0.32 × 1012$$\, \mathrm{M}_\odot$$, which agrees closely with other independent measures of the MW mass. We also discuss the MS in the context of H i streams in galaxy clusters, and find that the MS lies on the low-mass end of a continuum from Hickson groups to the Virgo cluster. As a tracer of the dynamical mass in the outer halo, the MS is a particularly valuable probe of the MW’s potential. 

Abstract We describe the discovery of a solar neighborhood ( d = 468 pc) binary system with a main-sequence sunlike star and a massive noninteracting black hole candidate. The spectral energy distribution of the visible star is described by a single stellar model. We derive stellar parameters from a high signal-to-noise Magellan/MIKE spectrum, classifying the star as a main-sequence star with T eff = 5972 K, log g = 4.54 , and M = 0.91 M ⊙ . The spectrum shows no indication of a second luminous component. To determine the spectroscopic orbit of the binary, we measured the radial velocities of this system with the Automated Planet Finder, Magellan, and Keck over four months. We show that the velocity data are consistent with the Gaia astrometric orbit and provide independent evidence for a massive dark companion. From a combined fit of our spectroscopic data and the astrometry, we derive a companion mass of 11.39 − 1.31 + 1.51 M ⊙ . We conclude that this binary system harbors a massive black hole on an eccentric ( e = 0.46 ± 0.02), 185.4 ± 0.1 day orbit. These conclusions are independent of El-Badry et al., who recently reported the discovery of the same system. A joint fit to all available data yields a comparable period solution but a lower companion mass of 9.32 − 0.21 + 0.22 M ⊙ . Radial velocity fits to all available data produce a unimodal solution for the period that is not possible with either data set alone. The combination of both data sets yields the most accurate orbit currently available. 

ABSTRACT We show that smoothed particle hydrodynamics (SPH) simulations of dwarf galaxies interacting with a Milky Way-like disc produce moving groups in the simulated stellar disc. We analyse three different simulations: one that includes dwarf galaxies that mimic the Large Magellanic Cloud, Small Magellanic Cloud, and the Sagittarius dwarf spheroidal; another with a dwarf galaxy that orbits nearly in the plane of the Milky Way disc; and a null case that does not include a dwarf galaxy interaction. We present a new algorithm to find large moving groups in the VR, Vϕ plane in an automated fashion that allows us to compare velocity substructure in different simulations, at different locations, and at different times. We find that there are significantly more moving groups formed in the interacting simulations than in the isolated simulation. A number of dwarf galaxies are known to orbit the Milky Way, with at least one known to have had a close pericentre approach. Our analysis of simulations here indicates that dwarf galaxies like those orbiting our Galaxy produce large moving groups in the disc. Our analysis also suggests that some of the moving groups in the Milky Way may have formed due to dynamical interactions with perturbing dwarf satellites. The groups identified in the simulations by our algorithm have similar properties to those found in the Milky Way, including similar fractions of the total stellar population included in the groups, as well as similar average velocities and velocity dispersions.